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Abstract 

The canonical front form Hamiltonian for non-Abclian SU(N) gauge theory in 3+1 dimen- 
sions is mapped on an effective Hamiltonian which acts only in the Fock space of one quark 
and one antiquark. The approach is non-perturbative and exact. It is based on Discretized 
Light-Cone Quantization and the Method of Iterated Resolvents. The method resums the dia- 
grams of perturbation theory to all orders in the coupling constant and is free of Tamm-Dancoff 
truncations in the Fock-space. Emphasis is put on dealing accurately with the many-body as- 
pects of gauge field theory. Pending future renormalization group analysis the running coupling 
is derived to all orders in the bare coupling constant. — The derived effective interaction has 
an amazingly simple structure and is gauge invariant and frame independent. It is solvable 
on a small computer like a work station. The many-body amplitudes can be retrieved self- 
consistently from these solutions, by quadratures without solving another eigenvalue problem. 
The structures found allow also for developing simple phenomenological models consistent with 
non-Abelian gauge field theory. 
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1 Introduction and Motivation 



One of the most important tasks in hadron physics is to calculate and understand the mass spec- 
trum and the wave functions of physical hadrons from a covariant quantum field theory. The wave 
functions encode all the properties which are needed for a phenomenological description of exper- 
iments. Wave functions and probability amplitudes have a natural appearance in a Hamiltonian 
approach. 

But the Hamiltonian bound state problem is notoriously difficult in a field theory. Procedures 
like those of Schwinger and Dyson or of Bethe and Salpeter are not easy to cope with in practice as 
reviewed recently in . Actually, the difficulties were considered so enormous that the Hamiltonian 
approach was given up in the Fifties altogether in favor of Feynman's action oriented approach. 



3| govern the scene, maturing from infancy 
are closer to experiment and have the 



Modern successors like the Lattice Gauge Theories 
in recent times §. Phenomenological models |5|, |6|77, 
different objective to classify the bulk of empirical data. They leave little doubt that a heavy meson 
contains primarily a pair of constituent quarks and not an infinity of sea particles as suggested by 
the perturbative treatment of quantum field theory. 

How can one reconcile these models particularly the constituent quark model (CQM) with the 
quantum field theory of chromodynamics (QCD)? There are several reasons why the front form || 
of Hamiltonian dynamics is one of the very few candidates, see [11, 13[ or [14|. Particularly 
the simple vacuum and the simple boost properties |l3|, 15, [l6|, 17 1 confront with the 'complicated 



vacuum' and the 'dynamical boosts' of the conventional approach, the instant form. These aspects 
are stressed also in the arguments of Wilson 18] and collaborators [[OJ. Nowadays we know that 
even the front form vacuum is not simpl e [p0|, pi] , 22 1, but still simplier than in the instant form: 
The problem can at least be formulated |23| , pMf . 

The success of discretized light-done quantization (DLCQ) particularly in 1+1 dimensions has 
stirred hope that its apparent simplicities carry over to the physical 3+1 dimensions j 



27, 



But this meets problems, among them: The Hamiltonian matrix increases exponentially fast with 
the particle number of the Fock states and with the number of transversal momenta. Truncating 
the Fock space to 2 particles like Tamm [ pif l and Dancoff invokes perturbation theory, violates 
gauge invariance and generates non-integrable singularities. One must resort to ad hoc procedures 
to make things working [26]. Truncating to 3 particles [27, 28], the numerical results are inconclusive 
as to test an onset of convergence. 

Even worse than that: Thus far, it is obscure how any Hamiltonian including the one on the 
light cone could be subjected to a renormalization group analysis. This and the possibly large 
coupling constant has motivated Wilson |]l9[| to give up the one-to-one connection with canonical 
field theory and to propose a radically new procedure in which the renormalization properties of 
front form operators play the crucial role. The problem is a very deep one and appears also in 
Quantum Electrodynamics (QED). In QED, the smallness of the coupling constant obscures the 
fact that the simple Coulomb potential between two point charges has not been derived from field 
theory thus far with other than perturbative methods. Although higher order effects are expected 
to be small, and are indeed so, see for example 32, 33 1 and p4l], no manifestly non-perturbative 
and closed analytical procedures are at hand. This is remarkable. In QCD the problems are only 
accentuated due to the larger coupling constant. 

By practical considerations one must reduce the matrix dimension of the Hamiltonian. But 
the reduction to comparatively small matrices is more than merely a technical issue. It becomes a 
matter of principle. When analyzing the source of the difficulty one realizes that the many-body 
aspects of a field theory have not been mastered thus far: Precisely those are treated in perturbation 
theory and precisely those are responsible for the large matrix dimensions. It is here where the 
present work supposedly contributes. The many-body aspects are kept in the form as they want 
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to appear, as resolvents, and these are not expanded as perturbative series. 

Our point is then: One should solve first the many-body aspects of the canonical Hamiltonian, 
which has not been done thus far, and then think on the problems of renormalization. This is not 
in conflict with [ 19 1 , but we emphasize a different aspects of the problem. One is faced then with 
problems similar to conventional many-body physics, problems as one meets them in the theory 
of atoms, nuclei or of solids. Remarkably, one can carry out this programme essentially without 
assumption, and still arrive at a solvable equation. As a matter of fact it is also simple. 

At the core of this work is a new method, the 'method of iterated resolvents' to be introduced in 
section f|. It allows to do 'perturbation theory to all orders' and 'perturbation theory in medium' 
without a smallness parameter. For to be specific, some of the ingredients of earlier work [13], 
particularly the Lagrangian and the DLCQ-Fock-space for QCD, are summarized shortly in sec- 
tion Section || summarizes the theory of effective interactions as known from the literature [35] 
and displays why the Tamm-Dancoff approach |29], |3^[ is bound to fail. The results of section"^ 
are then applied to QCD in section [|. At the end of many formal manipulations one takes the 
continuum limit. The effective interaction of a quark and an antiquark will then turn out as a 
sum of two terms which have an intuitively appealing interpretation: (1) The effective potential 
U, generated by the exchange of one effective gluon which is absorbed either by the same or by 
the other quark; and (2) The effective annihilation interaction U a , where the quark and antiquark 
annihilate into two effective gluons. The potential U is derived explicitly in section |6], and section 
summarizes the results including a broader discussion on their possible future application. 



2 The front form Hamiltonian for Quantum Chromodynamics 



Both in non-relativistic quantum mechanics and in field theory, the Hamiltonian operator prop- 
agates the system in time. In a covariant theory the concept of 'time' can be generalized, since 
the space-time parametrization is arbitrary. But following Dirac §, there are no more than three 
standard forms how to choose generalized time and the corresponding Hamiltonian dynamics: the 
'instant', the 'front' and the 'point' form. In this section we shall summarize [|l(], |ll|, [12], 13 1 the 
front form Hamiltonian for QCD, actually for SU(N), 

The time-like coordinate is chosen as x + = t + z, the space-like coordinates as x = (x,y,t — z) = 
(x±,x~). The four- vector of space-time is thus x^ = x±, x~~). In QCD the vector potentials 



(A^) cc ' are 3x3 matrices, and the Dirac spinors carry a color index c, i.e. ty c 
The Lagrangian density is given by the hermitian operator 



1, 



TrF^F^ + - ^7"^ - m F )^ + h.c. 



with c = 1,2, 3. 



(1) 



expressed in terms of the covariant derivative D M = <9 M + igA^ and the color-electromagnetic fields 
F^ u = d^A" — d u A^ + ig [A^, A^]. The four components of the energy- momentum vector 



1 f 1 

P v = - dx~ d 2 x ± ( 2TrF +K F/ + - fe 7 +D^ + h.c 

2 Jn V 2 L 



g +v c 



(2) 



i.e. P u = (P + , P±, P~), are strict constants of the motion. Its space-like components P± = (P , P 2 ) 
and P + do not depend on the interaction and in momentum representation are diagonal operators. 
The time-like component P~ = 2P + depends on the interaction and propagates the system in the 
light-cone time x + , i.e. i-^rl^) = an d therefore is the proper front form Hamiltonian [13]. 

The contraction of these four operators 



pvp„ = p+p- 



= H^c = H 



(3) 
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is Lorentz invariant and refered to somewhat improperly but conveniently as the 'light-cone Hamil- 



tonian i?i,c' [13], or shortly H. One seeks a representation in which H is diagonal, 

H\%) = Ei\Vi) . (4) 

Note that the eigenvalues Ei and matrix elements of H somewhat unusually carry the dimension 
of an invariant-mass-squared. 

Periodic boundary conditions on C can be realized by periodic boundary conditions on the 
vector potentials A M and anti-periodic boundary conditions on the spinor fields ^ a because C is 
bilinear in the latter. One expands these fields into plane wave states and satifies the boundary 
conditions by discretizing the momenta, hence Discretized Light-Cone Quantization (DLCQ), i.e. 



P- 



Jn, with n = 5, §, • • • ,00 for fermion fields, 
Jn, with n = 1, 2, . . . , 00 for boson fields, 



7T 

and p± = — n±, with n x ,n y = 0, dbl, ±2, . . . , ±00 for both . (5) 

This is done at the expense of introducing two length parameters, L and L±, which define a 
normalization volume Q = 2L(2L±) 2 . More explicitly, the free fields are expanded as Fourier sums 

1 1 



* a (x) = ^ ]T ^= ( WP, A)e~^ + d\v a (p, A)e^) , 

V ' ' q V P 

and A,(x) = ^J2^{ a MP^)e- tpx + a\el(p,X)e^) , (6) 



particularly for the two transversal vector potentials A 1 = A]_, (i = 1,2). Each particle is on its 
mass-shell p^p^ = m 2 . Its four- momentum is p^ = {p + ,p±,p~) with p~ = (m 2 +p 2 )/p + - Each 
particle state il q" is then characterized by six quantum numbers: 

q = (n,n x ,n y ,X,c,f) = (p + ,p ± , A, c, /) = (x,p x ,\,c,f) . (7) 

The first three, (n,n x ,n y ), specify the space-like momentum and A the helicity j or j. A quark 
is specified further by color c and flavor /. Gluons carry no flavour, and the color index c is 
substituted by the glue index a. The creation and destruction operators like a q and a q create and 
destroy single particle states q, respectively, and obey the usual (anti-) commutation relations like 

As an advantage of DLCQ, all quantum numbers are discrete. One deals thus only with simple 
(and dimensionless) Kronecker symbols. The spinors u a and v a and the transversal polarization 
vectors e± are the usual ones fl0| and defined in Jl3| . 

A thorough treatment should include the zero modes of the gauge fields particularly those of 



A + . Their importance had been demonstrated [23, ^J] for the vacuum sector. Here one deals with 



the particle sectors and their massive excitations. Explicit calculations with and without them 
p7] yield however the same results in the continuum limit. The global and the gauge zero modes 



f23 , |24J are therefore discarded in the sequel, and the usual light-cone gauge [^, 13] A + = is used. 
The light-cone Gauss equation, i.e. d^F^^ = gJ^, see [13] and [24], and the expansions in Eq. (j6j) 



complete the specification of all vector potentials A''. The space- like integrations in Eq.(|2D can 
be carried out explicitly, leading essentially to Kronecker deltas. One ends up with the light-cone 
energy-momenta [13| as operators acting in Fock space, i.e. P u = P v (a q , a q ,b q ,b q , d q , d q ). The 



various terms in the Hamiltonian are conveniently classified by a sum of four terms |13[, i.e. 

H = T + V + F + S . (9) 
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The kinetic energy T survives the limit of the coupling constant g going to zero. Since it is diagonal 
in Fock-space representation, its eigenvalue is the free invariant mass squared of the particular Fock 
state. The vertex interaction V is the relativistic interaction per se. It is linear in g and changes 
the particle number by 1 and only by 1. Matrix elements of V which change the particle number 
by 3 (as in the instant form) are strictly zero in DLCQ: The vacuum does not fluctuate. The 
instantaneous interactions F and S are gauge artefacts, are consequences of working in the light- 
cone gauge and proportional to g 2 . The seagull interaction S conserves the particle number. The 
fork interaction F changes the particle number only by 2. As illustrated below in Figure |l|, each 
block in the Hamiltonian matrix is therefore either the zero matrix, or has only seagull-, or only 
vertex-, or only fork-interactions, with very simple matrix elements tabulated in |l3}| . 

There is one and only one reference state which is annihilated by all destruction operators, 
namely the Fock-space vacuum \vac). Therefore, the Hilbert space for the single-particle creation 
and annihilation operators is the Fock space. It is the complete set of all possible states 

|$,) = Nt ...b\ N « ... 4 <at 2 . . . a\„\vac) , (10) 

subject to be eigenfunctions of the space-like momenta, with fixed eigenvalues P + and P±, i.e. 

0-7T -> 

P + = Y.Pt = -rK, and P ± (11) 

v L V 

The sums run over all particles in a Fock state. As consequence of discretization, the Fock states 
are denumerable and orthonormal: (<J>j|3>j) = 5ij. As usual, the momentum fraction carried by the 
particle is denoted by x = p + /P + , and the sum of all fractions is constrained to J2u x v = 1- Note 
that the Fock states can be made color neutral. Since P + has only positive eigenvalues and since 
each particle has a lowest possible value of p + , the number of particles in a Fock state is limited for 
any fixed value of the harmonic resolution K [15, [l(|. Next to the simple vacuum, this is another 



pecularity of DLCQ. 

For to enumerate all possible Fock states for a meson with a fixed harmonic resolution, one 
needs a classification scheme. A possible one is displayed in Figure |l[ With Fock states being 
color-singlets, the lowest possible value one can have for K is K = 1. This allows for one qq-p&ir, at 
the most. For K = 2, the Fock space contains in addition the g g- and the g<jg-sectors. For K = 4, 
one has at most 8 particles, namely 4 qg-pairs. In the figure all 13 Fock-space sectors possible for 
K up to 4 are denumerated by n = 1, . . . , 13. Note that the classification of these sectors does not 
change when K is increased: one just adds more complicated Fock space sectors to the figure. The 
number of sectors grows quadratically with K and has the value Nk = (K + 1)(K + 2)/2 — 2. 

In analogy to the figure, one can rewrite Eq.(^|) as a block matrix equation: 

N K 

(n\H\m) (m|* f ) = E t (n\^i) for all n = 1, 2, ... ,N K . (12) 

m=l 

The numbers n (m) denumerate the sectors. The problem is solved if one can find the sector wave 
functions (n|\I/j) for one or several eigenstates Note that each sector contains many individual 
Fock states with different values of x, p± and A. 

Actually, right from the beginning one could have chosen the conventional procedure. One 
could have written Eq.(|i~2|) in the continuum limit, replacing sums by integrals according to 
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Figure 1: The Hamiltonian matrix for a SU(N)-meson. The matrix elements are represented by 
the letters S, V, and F, corresponding to seagull, vertex, and fork-interactions, respectively. For 
better orientation, the diagonal blocs are marked by (•) and the zero matrices by (■). (In the 
preprint this table is replaced by a figure with the diagrams.) 
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But apart from introducing innumerable integration variables and separate summation indices like 
color, flavor and spin, one tends to loose oversight in a set of coupled integral equations as implied 



by Eq.(12). Since there is nothing wrong about thinking like Dirac in terms of matrices, we continue 
to do so because we find it more convenient. 

Each sector with a given parton number has arbitrarily many Fock states whose transversal 
momenta add up to a given value of P±. One needs a cut-off. Following Lepage and Brodsky 
|1C, 11, 12] the Fock space is regulated by the condition that the free invariant mass shall not 



exceed a limit, i.e. 



^\ T , < An + • ( 14 ) 

The sector dependent mass scales A n will not be needed below. They govern how much the particles 
can go off their equilibrium values p± = and x v = rn u /J2i/ rn v The number of individual Fock 
states in DLCQ is thus finite, and the Hamiltonian matrix will have a finite dimension. A finite 
matrix can be diagonalized numerically. 

This - in a nutshell - is the programme of Discretized Light-Cone Quantization J15[ |. One tends 
to conclude that the solution of the bound-state problem for a gauge field theory like QCD is 
solvable on a computer. 



3 The Effective Interaction and the Tamm-Dancoff Approach 

DLCQ applied to gauge theory faces a formidable matrix diagonalization problem. The problem 
is even worse than in non-relativistic many-body theory. Here are the steps one should take in 
principle: In a first step one sets up the Fock space as discussed in the previous section. In a 
second step one calculates the finite dimensional Hamiltonian matrix as illustrated in Figure [|. In 
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a third step one diagonalizes the matrix numerically. Here is the bottleneck of the method: Sooner 
or later, the matrix dimension exceeds imagination, and other than in 1+1 dimensions one has to 
develop new tools by matter of principle. 

Intuitively one aims at something like an effective interaction between a quark and an antiquark, 
similar to the effective interaction between a negative and a positive point charge. Effective interac- 
tions are a well known tool in many-body physics |3j|. To the community the method is known as 
the Tamm-Dancoff- Approach, as applied first to Yukawa theory for describing the nucleon-nucleon 



interaction [29, 30 1. A look on its salient features is worth the effort. 

Consider a Hamiltonian matrix as in Eq.(|l2"|) subject to diagonalization, H\^) = E\^). The 
matrix dimension be N. Explicitly written out, the eigenvalue equation reads 

N 

J2(i\H\j)(m=E <i|tf) , fori = 1,2,..., iV. (15) 

3=1 

The rows and columns of the matrix can always be split into two parts. One speaks of the P-space 
P = J2j=i \ with 1 < n < N, and of the rest, the Q-space Q = 1 — P. Eq.([l5|) can then be 
rewritten conveniently in terms of block matrices like 

(p\h\p) (p\h\q}\ f(p\m_ F f(p\m hH) 

(Q\H\P) (Q\H\Q)) \(Qm) * \(Q\*)J ' 1 j 
or explicitly 

(P\H\P) (P|tf) + (P\H\Q) (Q|tf) = E{P\V), (17) 

and (Q\H\P) (P|tt) + {Q\H\Q) (Q|tf) = E (Q|tf> . (18) 

Rewrite the second equation as 

(Q\E - H\Q) (Q|¥) = (Q\H\P) (P\V), (19) 

and observe that the quadratic matrix (Q\E — H\Q) could be inverted to express the Q-space 
wavefunction {Q\*&) in terms of the P-space wavefunction (P|\J r ). But here is the problem: The 
eigenvalue E is unknown at this point. One therefore solves first an other problem: One introduces 
the starting point energy to as a redundant parameter at disposal, and defines the resolvent of the 
Q-space matrix as the inverse of (Q\lo — H\Q), 

G ^ - m;hw) ■ (20) 

With Eq.(|l9| ) one defines then 

<Q|¥) = <Q|*(u;)) = G q (lo)(Q\H\P) (P\V) . (21) 

Inserting this into Eq.([l7|) yields an eigenvalue equation in the P-space 

(P\H eS (u)\P) (P\* k (u)) = E k {oj) (P|**(w)) , with (22) 

(P\H eS (u)\P) = (P\H\P) + (P\H\Q) G Q (u) (Q\H\P) . (23) 

The effective interaction in the P-space is thus well defined: It is the original matrix (P\H\P) plus 
a part where the system is scattered virtually into the Q-space, propagating there by impact of the 
true interaction, and finally scattered back into the P-space. Every value of uj defines a different 
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Hamiltonian and a different spectrum. Varying uj one generates a set of energy functions E^{oj). 
Whenever one finds a solution to the fixpoint equation [38, E59j 

E k {u) = u, (24) 

one has found one of the true eigenvalues and eigenfunctions of H , by construction. 

One should emphasize that one finds this way all N eigenvalues of H, see also [38], irrespective 



of how small one chooses the P-space. Even if one chooses a P-space with 'matrix' dimension n = 1, 
the single energy function Ei(uj) contains the infomation on all eigenvalues: The N — 1 poles of 
E\(oS) generate the N eigenvalues by means of Eq.(24). Explicit examples for that can be found in 
|38| , |39f and in Appendix [B|. It looks as if one has mapped a difficult problem, the diagonalization 
of a huge matrix onto a simplier problem, the diagonalization of a much smaller matrix. This is 
true, but not completely. One has to invert a matrix. In practice, the inversion of a matrix is as 
diffficult as its diagonalization. In addition, one has to vary uo and solve the fixpoint equation (|24"|). 
The numerical work is thus rather larger than smaller as compared to a direct diagonalization. 

The advantage of working with the so defined effective interaction is that resolvents can be 
approximated systematically. The two resolvents 

Gq{uj)= (Q\.-T-U\Q) ' andg °H= { Q\ U - T \Q) ' (25) 

defined once with and once without the non-diagonal interaction in the Hamiltonian H = T + U , 
respectively, are identically related by 

G q {uj) = G (uj) + G (u)) U Gq(u) , or by (26) 
G Q (u) = G (u) + G (u)UG (uj) + G (uj)UG (lo)UG (uj) + . . . , (27) 

that is, by the infinite series of perturbation theory. The point is, of course, that the kinetic energy 
T is a diagonal matrix which can be trivially inverted to get Gq(u>). 

After these formal considerations we return to gauge theory and its Fock space representations 
as displayed in Figure [l]. The identification of the gg-space with the P-space and the related effective 
interaction between one constituent (q) with an other (</) appears as the most natural thing to do. 
But now the trouble starts: Like Tamm and Dancoff [p9| , [30| l one truncates the above series to 
the very first term H cS = H + H\Q)(Q\(u - T)- 1 |Q)(Q|F. But the Tamm-Dancoff procedure 
fails for light-cone kinematics: When taken literally, the effective interaction has a non-integrable 
singularity [26|. Even worse, simple estimates show that the series in Eq.(^Tj) diverges [27] order by 
order, all that irrespective of violating badly gauge invariance due to the truncation. One has to 
invent ad hoc measures to get things working p^]. One could try to resum the series in Eq.(p7|) at 
least partially, guided for example, like in Quantum Electrodynamics, by a smallness assumption on 
the coupling constant. When one tries that, one gets lost soon in the murky depths of perturbation 
theory driven to all orders, and quits the game. One the other hand, the procedure displayed in 



Eqs.fliq) to (26) is in principle exact. Maybe one should look at it from an other point of view. One 
of the many conclusions might be that the interaction term in the denominator of Eq.(p5|) should 
be dealt with in a better way, for example, by the 'method of iterated resolvents' to be presented 
next. 

4 The Method of Iterated Resolvents 

DLCQ applied to gauge field theory is, as mentioned, particular among matrix diagonalization 
problems to the extent that a finite number of Fock-space sectors like in Figure [j] appears in the 
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most natural way. The same holds true for quantum-mechanical ^4-body problems where the Fock 
space can be classified in terms of the finite number of O-particle-0-hole (0-ph), 1-ph, A-ph 
sectors, see for example [|40|[. In either case the block matrix in sparse. Most of the blocks are zero 
matrices due to the nature of the underlying Hamiltonian. One should make use of that! 

In DLCQ one is confronted with the diagonalization of a finite dimensional block matrix, for 



which Eq.(16) represents the example of 2 x 2 blocks. The step from Eq.(16) to Eq.(|22|) can then 
be interpreted as the reduction of a block matrix dimension from 2 to 1. As a matter of fact, if 
one chooses the Q-space identical with last sector Nk in Figure and the P -space with the rest, 
Q = 1 — P. one reduces the block matrix equation from dimension Nk to Nk — lj with an effective 
interaction acting in the now smaller space. This procedure can be repeated, and again iterated, 
until one arrives at an effective interaction of block matrix dimension n = 1, with a well defined 
effective interaction in the 'l'-space. All that is needed is a reasonable notation [41]. 



Suppose, in the course of this reduction, one has arrived at block matrix dimension n, with 
1 < n < Nk- Denote the corresponding effective interaction H n (u). The eigenvalue problem 
corresponding to Eq,(|22|) reads now 



J2{i\H n (u)\j){j\^(u;)) =E{u) (i|tt(w)> , for* = l,2,...,n. (28) 

Observe that i and j refer here to sector numbers. Now, in analogy to Eqs.([20|) and (^3), define 

G n (uo) = 1 , and (29) 

{n\uj — H \n) 

n-l 

(n|*H) = G n (uj)^2(n\H n (uj)\j) 01*0*0) , (30) 

3=1 

respectively. The effective interaction in the (n — l)-space becomes then [41| 

H n ^(uj) = H n (u) + H n {uj)G n (uj)H n (uj) (31) 

for every block matrix element (i\H n -i(uj)\j) . To get the corresponding eigenvalue equation one 
substitutes n by n — 1 everywhere in Eq.(|28[). Everything proceeds like in section ||, including the 
fixed point equation E(u>) = uj. But one has achieved much more: Eq.([H]) is a recursion relation 
which holds for all 1 < n < iV#-! The only convention one must have, of course, is that one has 
started from the bare interaction H, thus Hn k = H. The rest is algebra and interpretation. 

How can the procedure be interpreted? One gets some inspiration from the paradigm of a 4 x 4 



block matrix as explicitly dealt with in Appendix ||. The consequences of Eq.(31) are worked out 
in Appendix [C]: When expressed in terms of the bare interaction (matrices) H and the resolvents 
G n , the effective interaction in sector n wants to develope chains like 

H n (u) = ... + HGi{u)HG m {u)HG r [uj)H + ... , with l,m,r > n. (32) 

The order of the chains is given by the number of propagators and not, as usual, by the power 
of the coupling constant. The maximum order of the chains is strictly finite for a finite matrix, 
as opposed to the infinite series of perturbation theory. The sector numbers in the chains are not 
arbitrary, but obey the rules tabulated in Appendix |^. For example, the effective interaction H n 
has no propagator with sector number n. This implies that the system can not fall back, a feature 
distincly different from the usual perturbative propagators displayed in Eq.(p5|). Moreover, most of 
the chains vanish from the outset for the DLCQ-Hamiltonian: It suffices that only one of the four 
bare block matrix elements appearing in Eq.(|32"|), namely (n\H\l), (l\H\m), (m\H\r), or (r\H\n) is 
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Figure 2: The Hamiltonian matrix for a SU(N)-meson when only vertex diagrams (V ) are included. 
Compare also with Figure [Ij. Zero matrices are marked by (■). (In the preprint this table is replaced 
by a figure with the diagrams.) 



K N p 


N p 

Sector n 


2 
1 


2 3 4 
2 3 4 


3 4 5 6 
5 6 7 8 


4 5 6 7 8 
9 10 11 12 13 


1 2 


qq 1 


• 


V 






2 2 
2 3 
2 4 


55 2 

qq a 3 

gggg 4 


V 


• V 

V • V 
V • 


V • • 

V 

V 




3 3 
3 4 
3 5 
3 6 


9 9 9 5 
qqqq g 7 

qqqqqq 8 




V 

V 

V 


• V • • 
V • V • 

• V • V 

• • v • 


V ... 

V 

V 

V 


4 4 
4 5 
4 6 
4 7 
4 8 


9 9 9 9 9 

qq g g g io 

qqqqqq g 12 
qqqq qqqq 13 






V • • 

V 

V 

• • V 


• V 

V • V • 

• V • V ■ 

• V • V 
V • 



a zero-matrix. Unpleasant is to some extent, that one has to deal with 'resolvents of resolvents', 
or with 'iterated resolvents'. On the other hand, each resolvents represents a resummation of all 
orders in perturbation theory, and it is this particular aspect which makes them interesting, at 
least. 



5 Application to Quantum Chromodynamics 

The method of iterated resolvents is almost ideally suited for dealing with the DLCQ matrix for 
Quantum Chromodynamics, see Figure |lj, and for deriving the effective Hamiltonian in the gg-space. 
In doing so, we actually shall use a trick which will simplify considerations enormously. Practi- 



tioners in Light-Cone Time-Ordered Perturbation Theory |10| , 11, 13, 25] know that they can omit 
the instantaneous interactions F and S until they actually compute a particular diagram. Then, 
every intrinsic line in a graph must be combined with the instantaneous partner line associated 
with the gauge artefacts. Only then, the sum of all time ordered diagrams becomes manifestly 
identical with gauge-invariant Feynmann scattering amplitudes. There is no exception known to 



this rule, thus far, in all graphs computed explicitly. In the sequel, this 'gauge trick' [25] is adopted 
by setting formally to zero all block matrices in Figure [l] which are gauge remnants. One then gets 
a block matrix structure as displayed in Figure ||. The extreme sparseness of this block matrix will 
allow us to carry the procedures through to the end. When eventually arriving there, we will check 
explicitely in Appendix [D] whether the restitution of the gauge artefacts changes anything in the 
structure of the solution. The answer will be: 'It doesn't'. 

When writing down the chains for the block matrix displayed in Figure ^ according to the rules 
of the method of iterated resolvents, the effective Hamiltonian in the gg-space becomes 

H x = Ti + VG 3 V + VG 3 VG 2 VG 3 V . (33) 
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1 1 



3 

Ui2 



Figure 3: Three graphs of the effective inter- 
action in the qq-space. — The lower two graphs 
correspond to the chain VG3V, the upper cor- 
responds to VG3VG2VG3V . Propagators are 
represented by vertical dashed lines with a sub- 
script 'n ' for the sector. 





Figure 4: Graphs of the spectator interac- 
tion in the qq g-space. Note the role of the 
gluon as a spectator. 



Only two chains survive the procedure, one with one, and one with three propagators Gj = Gj(u). 
One checks explicitly with Eqs. flll4| ) and (115) that no chains with four, five or more propagators, 
are possible, even not for the case of arbitrarily large harmonic resolution K. Eq.(|33"|) stands thus 
for full QCD in the continuum limit. The gauge trick is thus a great help: Including the gauge 
artefacts, the effective Hamiltonian in the qq-spa.ce would given 32 non-trivial chains up to order 
4. These are explicitly tabulated in Appendix ^. 

Why are there only two chains in Eq.(p3|)? The system can scatter out of the gg-space into the 
qq g-space by emitting a gluon from the quark or from the antiquark. There is no other way since 
all other block matrix elements (l\H\j) vanish according to Figure ||. This kills all possible chains 
which do not start like VG^VGjVGi .... The scattering from the qq g-spa,ce proceeds either by 
scattering back to the qq-space, which is the first chain in Eq.(33), or by scattering into another 
space by a sequence like G^VGj. The rules of MIR require 3 > j > 1. The only solution j = 2 
gives the second chain. 

What do the chains in Eq.(]33|) mean in terms of physics? We find it convenient to illustrate 
that by the diagrams in Figure ||. The first term in the effective Hamiltonian Hi refering to the 
kinetic energy (7\) in the qq-space needs no illustration. In the second term (VG3V) the vertex 
interaction V creates a gluon and scatters the system from the qg(l)-sector into the qq g-space. 
As indicated by the vertical line with the subscript '3', the three particles propagate there in 
'3'-space under impact of the full interaction before the gluon is annihilated. The gluon can be 
absorbed either by the antiquark or by the quark, corresponding to the two 'interactions' Ui t i and 
U1.2, respectively, in the figure. Two further graphs, describing the irradiation of the gluon by 
the antiquark, are not shown; like throughout in the sequel we use the graphs for the purpose of 
illustration rather then of a complete description. As the inverse of the non-diagonal Hamiltonian 
#3(0;), the propagator G^{u) is non-diagonal as well, in general. Therefore, despite the total (space- 
like) momenta being conserved strictly by the interaction and thus by propagation, the space-like 
momenta of the individual particles can change. A finite size box instead of a thin line would 
thus be more appropriate for graphically representing propagators. With this proviso in mind, and 



12 




Figure 6: Twelve graphs corresponding 
to Fig. with the gluons absorbed after 
propagation. 



the recognition that an interaction can have several graphs, we proceed. The last term, finally, in 
Eq.(|33"|) is represented correspondingly by the graph U\fl in Figure ||. 

5.1 The quark-pair-glue propagators, and propagation 'in medium' 

The mind trained by perturbation theory wonders about all those zillions of graphs, he is used to 
work with. The answer is that they have neither been ommitted nor neclected but that they reside 
in the propagators Gi, very effectively 'resumed to all orders'. How do they look? 

By the rules established in Appendix |C], one obtains consecutively for the effective Hamiltonians 
in the sectors with one qq-pair and 1,2,3 or more gluons 

H 3 = T 3 + VGqV + VG 6 VG 5 VG 6 V + VG 4 V , (34) 
H 6 = T 6 + VG W V + VG W VG 9 VG W V + VG 7 V , (35) 
H w = Tjjo + VG l5 V + VGisVGuVGinV + VG n V , (36) 

respectively, to whose resolvents we refer collectively as the quark-pair-glue propagators. Note that 
these are all exact relations, by the same reasons as discussed above, for the qq(l)-sector . Note also 
that they all have the same structure. They differ from the effective Hamiltonian in the qq-sp&ce, 
Eq.(|33"|), only by their respective last term, in which a gluon is annihilated into a qq-p&ir. Since the 
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gg-sector has no gluon, by definition, the corresponding term must be absent there. In passing, one 
notes the 'replica-structure' of these Hamiltonians and propagators. Like in Russian puppets one 
has propagators of propagators, and a certain amount of self-similarity cannot be denied. These 
strucures seem to be particular to gauge theory, abelian or not. 

Let us have a look at the diagrammatic representation of these Hamiltonians, and in particular 
on the effective Hamiltonian in the qq g(3)-sector, Eq.(34). A typical collection is given in Figures |] 



and [j| The subset of graphs displayed in Figure |] looks exactly like the effective interaction in 
Figure ||, except of the additional gluon. The gluon does not change quantum numbers under 
impact of the interaction, and acts like a spectator. We thus refer to the graphs in Figure |4] as 
the 'spectator interaction' U3. In the graphs of Figure ||, the gluon undergoes scatterings which 
correspondingly are refered to as the 'participant interaction' U3, in the sequel. The effective 
Hamiltonians the quark-pair-glue sectors qqg . . . g, as given in Eqs.(|34|) to (|36|), can thus be written 



H n = T n + U n + U n , for n = 3, 6, 10, 15, . . . , (37) 

collectively, and with the obvious corresponding definitions. 

In terms of physics, the separation into spectators and participants becomes more transparent 
in a perturbative consideration. Think of the incoming gluon in Figure [| as being irradiated by a 
quark, and of the outgoing gluon as to be absorbed either by the quark (a) or by the antiquark 
(b), and draw the corresponding diagrams. Each diagram in Figure ||| will then be associated 
with two diagrams in Figure |(| The labels on the diagrams in the figures have no other purpose 
in the present context than to help the reader with this association. The diagrams in Figure ^ 
look familiar. Would one do plain perturbation theory, they would be diagrams of order 4 in the 
coupling constant. Indeed, they look like the familiar mass, vertex or vacuum 'corrections', which 
are associated usually with the running coupling constant |45|, Q. Here, however, we continue to 
deal with propagators in medium and take the analogue only as an argument why the separation 
into 'spectators' and 'participants' might be useful. 

The essence is that the propagators G n associated with the spectator interaction U n , i.e. 

G n = — , ( while G n = — =- = l — \ , (38) 

U-T n -U n ' \ U-T n -U n -U n U ~ H n J ' 17 

are identically related to the full propagators G n by 

G n = G n + G n U n G n . (39) 
As usual, they can be written as an infinite series 

G n = G n + G n U n G n + G n U n G n U n G n + . . . . (40) 



The main difference to the usual series like in Eq.(p6[) is, that there the 'unperturbed propagator' 
Gq(u) refers to the system without interactions while here the 'unperturbed propagators' G n contain 
the interaction in the well defined form of U n . One therefore deals here with 'perturbation theory 
in medium'. Note that the present series is different from the above Eq.(p6[) also with respect to 
the physics: The system stays in sector n. This allows for an identical rearrangement of the series 



G, 



I _ 3 _ _ 

1 + ~zG n U n + —G n U n G n U n + . . . 
z o 



G r . 



1 _ 3_ _ 

1 + ~zU n G n + — U n G n U n G n + . . . 

2 o 



(41) 



which can be verified order by order and which, to our recollection, has not been given before. The 
series in the square bracket are known to be the inverse square, i.e. 

Rn = ! l _ = =l + ]p n Un + fe n LZ„ G n U n + . . . , (42) 
1 — G n U n 
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and therefore the full propagators can be rewritten identically as 



Gn — RnG n Rn — . = G n . = . (43) 



l-U n G n ^l-G n U n 
This result can be verified also by closed operator relations. First, one notes the strict identity 

(1 - U n V n )(u> - H n ) = (l- U n 1=-) (u-H n )=u- H n . (44) 



'u - H n 



Since its adjoint is also correct, one gets 



GlG n = 1 —— ■= i . (45) 

1 - U n G n {uj - H n f 1 - G n U n V ' 



This is a positive valued operator. Taking the square root with this precaution, gives Eq.(|43|) as an 
identity. One thus can proceed without having to worry on the square roots, their ambiguities of 
signs, and on the convergence of the series. The square matrix R will always be sandwiched between 
a quark-pair-glue propagator G and two vertex interactions V, for which reason it is convenient to 
introduce V as an abbreviation, defined by 

VG n (u) V = V Ri(u)G n {u)R n {uj) V = VG n V . (46) 

Below, a very natural and physical interpretation is given to the operator R, as being related to 
'running coupling constant', but here we continue to proceed formally. We use the above findings 
to systematically rewrite Eqs.(|33|)-(|3"6|). Working upwards in the hierarchy, one gets consecutively: 



H 6 = T 6 + VG W V + VG W VG 9 VG W V , (47) 
H 3 = Ts + VGeV + VGeVG^GeV, (48) 
H 1 = H 1 = Ti + VG 3 V + VG 3 VG 2 VG 3 V . (49) 

Instead of being similar, the sector Hamiltonians H n = T n + U n become equal for suffiently large 
K. All what is different is that they act in different quark-pair-glue spaces, namely in the gg-spaces 
with 1, 2 or more gluons. 

The gluons, however, by construction do not take part in the interaction. They act like inert 
spectators, very much like displayed already in Figures | and | for the interactions in the qq(l)- 
and the qqg-space, respectively. All what one has to do is to update the figures and to replace each 
point-like vertex by say a square, which should symbolize the impact of the operator R. 



5.2 The key point: The quark-pair-glue resolvents 'in the solution' 

The operators H n = H n {uo) are bona fide Hamiltonians which can be diagonalized on their own 

H n (u)\* n (u)) = , for n= 1,3,6,10,... . (50) 

Since the (spectator) gluons do not interact with the bound qq system, by construction, one can 
relate the various spectra, knowing nothing on U n except that the same potential acts only between 
the quarks. For n = 1, the Fock states are \q; q) = b^d^vac) and the eigenvalue equation reads 

Y J {q)q\Hi{^)W^')W^\^)) = MZ(u)(q;q\* b (u;)) . (51) 

q',q' 

The eigenstates labeled by b form a complete set of states and denumerate the physical mesons 
and their spectrum. They correspond to physical particles up to the fact that their invariant mass 
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Mb = Mf,(u) still depends on the starting point energy. Correspondingly for n = 3, one has 
\q;q;g) = b\d\a\\vac) and 

]T (q;q; g\H 3 (u)\q'; <f\<f)W\<f\cf\1> b ,s{u)) = Ml s {u){q;q;g\^ btS {u)) . (52) 
q',q',9' 

Because the gluon is a non-interacting spectator, the eigenfunction is a product state |i/>j, s (o;)) = 
\^b(u})) <g> \<p s ), labeled by the two indices b and s. In fact, the eigenvalues in this sector must be 

Ml = + % , (53) 

V 9! 9 

corresponding to the free gluon moving with momentum (fraction) k g± (x g ) relative to the meson 
with mass and momentum (fraction) kb ± = —k g± (x& = 1 — x g ). Correspondingly the eigenvalue 
in the qqg . . . g sectors are 



M 2 + k 2 a 2 +k 2 
i-i _ m b + K g± _j_ /V-g + K g± 

(1 x^) 



M^g g = 1 _ ^ + ^••• g _ i g " . (54) 



Here the meson moves against the cluster of free gluons with free invariant mass f-i g ...g- Note that 
these statements are frame independent and a direct consequence of the front form. Only in the 
front form the transition to a moving frame is trivial. In the the instant form the boost operators are 
non-diagonal and complicated. Note also that these statements hold if and only if the interaction 
in all sectors is identical, i.e. 

U n {u) = Ux{u) . (55) 

Is this really true? A word of caution seems to be in place. Actually, it is not true for any finite 
K. The argument is simple: The 'continued fraction' behind the resolvents of resolvents and the 
Russian puppet structure of the theory becomes interrupted in the last sector. For a finite K there 
is always a last one, no matter how large K is chosen, and thus the continued fractions in the 
different sectors do not have the same length and therefore are different. In the continuum limit, 
however, this argument of denumeration does not hold. Thus, for any finite but sufficiently large 
K the above statement is only 'sufficiently true'. 

The resolvents of H n , in general, are non-diagonal matrices in Fock space representation. Trans- 
forming to the representation in which H n is diagonal, one gets for n = 3 

(q;q;g\G3(u>)\q';q;g') = {q;q;g\ , . \q';ef;g f ) 

uj - H 3 {u) 

= ^(^sII'mX^mI — wr • (56) 

b,s uj - b g f - -EL 

(l — Xg) Xg 

Consider now the resolvent G(cj) as function of u. It depends, of course, on the value of u and 
on the spectral distribution of the eigenvalues. Henceforward, as being reasonable for a 'bound 
state calculation', we shall restrict to the spectral region where the eigenvalues are sufficiently well 
separated from each other. G{oS) will be strongly peaked whenever u hits one of the true eigenvalues, 
say u = M 2 ,(oj), according to Eq . (p4|) . At the peak of this function, however, the following sequence 
of steps can be performed: 

(q;q;g\G 3 (uj)\q';q';g') = ^fag;^,,) TTTTti fT (^b,s\q'\ 4\ d) (57) 

bs M 2 b g ^ 91 - 

M b> ~ 71^7 1- 



16 



1>; q; g\* b ,s) ^-(^W; 4\ a') (58) 

6,* Ml - Ml- 

O Xg 

X 



5 -E<95 9^I*mX*mI9 / ;9';^> (59) 



A; 2 



= -g- (q;q;g\q';q';g') . (60) 

The first step invokes a new; smallness condition: 

x g <^e<l and fc| ± < A 2 < M fe 2 . (61) 

It has a completely different origin than the cut-offs A n introduced for Fock-space regularization. 
The second step, including closure, is not only permitted but also a valid approximation in line 
with standard many-body technology. 

No approximation, however, is involved relating the off-shell mass k 2 ± /x g of the gluon to the 
other particles in \q;q;g). Denoting the single-particle four-momentum of the gluon, the quark, 
and the antiquark by k g ,k q , and kg, respectively, only well-known light-cone kinematics is needed 
to establish that the two four-vectors 

^=(k g + K-^T = ^^ and F q = (k g + k q -k' q Y = (62) 

are identical, and proportional to the time-like vector 

< = ( ?? + ,r/ ± ,r ? -) = (0,0 ± ,2) , (63) 

with rj 2 = 0. The identities 

^ = (k g + k q + k' q } 2 -(k q + k q ) 2 = - — (k'^-ka) 2 (64) 

Xg y Xg V y / 

= {kg+k q +k' q ) 2 -(k' q + k' q ) 2 = -— (k q -k' q ) 2 (65) 

Xg 

are then obtained straightforwardly. 

Some words of caution seem to be in place. (1) The smallness condition requires that A is 
smaller than the lowest bound state. Restraining the transverse momentum of the gluon to be less 
than the pion mass is not a serious draw-back. One should note that the off-shell mass of the gluon 
can be large irrespective of how well the smallness condition is fulfilled. (2) One can state with 
certainty, that the above reduction is inadequate if not false at sufficiently large values of u, where 



one might be in the region of 'overlapping resonances' [42]. This will require different technologies 
and possibly those based on random matrix models p3| , 44]. It is there where the concept of a 
temperature will find its way into the bound state problem of gauge field theory. 

Finally, the resolvent in the qq ^-sector has determined itself to be very simple: In the solution, 
the resovents G n are diagonal in Fock-space representation and strictly independent of u. Instead 
of having 'resolvents of resolvents' the hierarchy of resolvents is broken. The particles 'in medium' 
propagate like free particles. 



6 The effective interaction for Gauge Theory 

The above considerations are summarized as follows. The effective Hamiltonian 

H cS = T + U + U a (66) 
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Figure 7: The ful/ effective interaction in the 
qq-space, resumed to all orders in the bare cou- 
pling constant. Circles represent the running 
coupling constant. The one-gluon exchange in- 
teraction U = Ui t i + U\ : 2 provides the binding. 
The annihilation interaction U a = can be 
active only if the quarks have the same flavor. 
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Figure 8: Two typical graphs, taken from Fig- 
ure ||, should illustrate how the propagator box 
(dashed line) contracts into a square around 
the vertex. 



is supposed to be diagonalized only in the qq-space. The eigenvalue equation H e g \ipb) = E |^&) 
can be understood as the finite dimensional (N s ) matrix equation 

N s 

J^(i\H ee \j)(j\ip b ) = E b (# 6 ) , fori = l,2,...,N s . (67) 

Subject to fixed values of P + and P±, the basis states |j) denumerate all possible Fock states 

\j) = = &g4l )' or \<li0 = \ x >k±;X q ,\q) ■ ( 68 ) 

It often suffices to label an individual Fock state by the momenta of the quark x = x q and k± = , 
since the antiquark has Xq = 1 — x and = —k±. The effective Hamiltonian has a kinetic energy 
T and two kinds of potential energies U. As operators acting in Fock space they are defined by 

T = £^=f±^ + £^p±£) f , (69) 

u = H b\d\dq,b ql {q,q\U\i/,tf) 5 fqJql Sf 9 j_, , (70) 
q,q,q',q' 

U a = b 54 d ?'V U a \q,q) S fqtfq 5 fq ,j_, . (71) 

q,q,q',q' 

One distinguishes a flavor-conserving (U) and a flavor-changing effective interaction (U a ). Both 
of them scatter a quark from state q to state q', and the antiquark from q to q' . They are given 
diagrammatically in Figure [?], and more formally by 

U = VG 3 V = VRlG 3 R 3 V , and (72) 

U a = VG 3 VG 2 VG 3 V = VR\G 3 R 3 VG 2 VRIg 3 R 3 V . (73) 
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Most importantly, the propagator G 3 (but not G 3 l) 'in the solution' is diagonal in Fock space 
representation, its numerical value being closely related to the momentum transfer Q of the quarks, 

(q;q;g\G 3 \q";q";g") = ^ {q;q;g\q";q";g") , with Q 2 = (k q - k' q f . (74) 

The operators R had been defined in Eq.(|42|); here they read 

R 3 = L = . (75) 

V 1 - G3U3 

What is their interpretation? 

6.1 The Running Coupling Constant 

The role of R becomes more obvious when setting R = 1 but keeping R^ in Eq.(fT^), 

U ~ VRlG 3 V . (76) 

Expanding R) to first non-trivial order, picking out a particular term, one gets diagrams like in 
Figure ||, or formally 

U ~ VG 3 V + ^VG 3 VG 6 VG 3 V + . . . . (77) 

As mentioned already, a propagator should be represented rather by a box than by a vertical 
line. This becomes important when the propagator is diagonal, see Eq.(74|), i.e. when it does 
not change the momenta of the 'in-' and 'out-states'. As shown in Figure 8[ the surrounding box 
lines can then be contracted over the exchanged gluon. The in- and out-momenta of the R-hox. 
are then determined only by the quark single-particle four-momenta k q and k' q . Note that these 
two also determine uniquely the four-momentum of the irradiated gluon. Therefore, in Fock-space 
representation, the operator R must be a function of k q and k' q alone, and thus only a function of 
the four-momentum transfer Q as defined in Eq. ([74]) . The vertex function 



(k q \R\k' q , k g ) = r(k q , k' q ;A, m q ,g) = r(Q, A) (78) 

depends thus on all variables in the problem, particularly on Q and the cut-offs A. As the figure 
demonstrates, the operator R and its matrix elements are closely related to the running coupling 
constant: both have the same perturbative expansion |^, E^]. But it should be emphasized that 
perturbative expansions are used here only for the purpose of illustration. The present formalism, 



particularly Eq.(75), resumes it to all orders in the bare coupling constant, contrary to the familiar 



expressions gained from perturbation theory 46]. One should mention here the potential danger 
hidden in the cut-off depence of r(Q, A). Perturbative estimates yield diverging expressions for ever 
increasing value of A. The problems related to these well-known divergencies will be discussed in 
future work, see also section fO- 



6.2 The effective interaction 

The actual computation of the vertex function r(Q,A) as function of all its parameters clearly goes 
beyond the scope of the present work and must be left to the future, despite its importance. By 
the same reason, the computation of the annihilation interaction can not be given here. Rather 
shall we restrict ourselves here to calculate explicitly the flavor-conserving interaction U. 
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Consider the matrix element U\j_ in Figure ^. In this graph, an initial state \q,q) state is 
scattered into a final state \q',q'), going through an intermeditate state \q', q, g). As an adventage 
of DLCQ, one can introduce these states as to be ortho-normalized and invariant under SU(N): 



\q,q) 



\q ,q,9) 



Y^bl{k q ,\ q )d\{kq,Xq 



c=l 



\vac) 



\vac) 



and 



c=l 



12 T c,c'bUk' q ,X q )dt(k q ,Xq)al(kg,X s 



vac) 



(79) 
(80) 

(81) 



c=lc'=l a=l 



The color matrices T" c , are normalized as usual by 2Tr(T a T ) = 8^. Since the gluon has a positive 
value of kg, one necessarily has either 



or 



(82) 



We restrict here to the former. The calculation of the vertex interaction proceeds in two steps. 
By means of the tables [13], one calculates the vertex matrix elements, separately for emission and 
absorbtion of the gluon, 



(q,q\ VR ] \q',q,g) 
W,q,g\RVW,<?) 





/<?2 r 2(Q,A) 


/ 2n r V 





u(k q , X q ) j^e^kg, X g ) u(k' q , X' q ) 



/ n 2_! / 9 2 r 2(Q,A) 



u(k q , X q ) j"e*(k g , Xg) u(k' q , X' q ) 



2n r 



(83) 
(84) 



respectively. Since both the (LC)-Hamiltonian and the Fock states are color neutral, all color 
algebra has reduced into the overall factor (r? c — l)/2n c . According to Eq.(72) one has to sum over 
all intermediate states. Since the momenta are fixed, only the sum over the two gluon helicities 
remains, which is done in the usual way: 



(85) 



The time-like vector rj had been defined in Eq.(62). In Appendix |E| it is shown how the gauge 
remnant of the polarization sum cancels exactly against the contribution from the instantaneous 
interaction. The latter, remember, appears since one works in the light-cone gauge A + = 0. The 
remainder due to g^ u becomes then manifestly gauge invariant: 



{q,q\u\q',q' 



2{2itf 



{q,q\u\q' ,q') 



, - lrrl , 1 /9(0,A) (X q ,X q \S(Q)\Ki X '«) 

with {q,q\U\q,q) = -— =■ 



4vr 2 Q 2 



\Jxq{\ X q ) X q (l X q ) 



(86) 
(87) 



and {Xq,X q \S(Q)\X q ,X'q) 



u(kq, X q ) 7 M U(k q , X q ) V,(k q , X q ) 7^ u(k q , X' q 



The spinor factor S(Q) collects the familiar and manifestly Lorentz invariant current-current in- 
teraction, which ultimately will be responsible for the fine and hyperfine structure. 
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The bare coupling constant g is combined with the vertex function r{Q, A) into the like-to-be 
'running coupling constant' (3(Q,A), 



/?(Q,A) 



n,. 



1 9 2 



2n r Airhc 



3 (Q,A) 



For QCD the gauge group factor {n 2 — l)/2n c reduces to 4/3, for QED it has the value 1. It is 
here where the essential difference of Abelian and non-Abelian gauge theory appears: The coupling 
constant runs differently. As a consequence, the one is confining and the other one is not. — The 
graph Ui t 2 in Figure [?] is computed correspondingly. 

In the discretized case of the matrix equation, no divergencies can occur because of the reg- 
ularization. In practical work, however, and by matter of principle, it is more convenient to go 
to the continuum limit K — > oo, simply by replacing sums with integrals. The conversion factor 
from Eq.(|i~3"l) cancels against the corresponding factor in Eq.([S6l), and the matrix equation (|67|) is 
converted into the integral equation 



M^(x,k±;Xq,Xg\ip b ) 



m a + kf m g + kf 



+ 



x 



(x,k±;X q ,Xq\ip b ) 



+ 4tt 2 ^2- 



Q 2 



\J x gO- x q) x 'qO- x 'q) 



(x, k±; X q , Xq\ip b ) ~ ( x ', k±; X' X'M b ) . (90) 



The A should remind to the fact that the integration over the perpendicular momenta is restricted 
to a range consistent with Eqs.(|i~4l) and (|6lt). Up to the omitted annihilation term, the resulting 
eigenvalues M 2 must coincide with the eigenvalues of the full Hamiltonian as shown above. The 
mass Mb must be interpreted as the mass of a physical meson. The corresponding wavefunctions 
{x, k± ; X q , Xq\ipb) gi ve the probability amplitudes for finding in that meson a quark with momen- 
tum (fraction) k± (x) and helicity X q . Finally, one should emphasize that the front-form effective 
Hamiltonian, appearing as the kernel of the integral equation (p0|), is manifestly boost invariant. 



6.3 Renormalization Group analysis 



The integral equation (|90| ) can be solved for any fixed value of the cut-offs A. Therefore, both the 
wavefunctions and the eigenvalues depend explicitly on it, Mb = Mb(A). 

This is unphysical, since a physical results should not depend on mathematical tricks like 
regularization. Beyond that it is potentially dangerous since the vertex function diverges as a 
function of A. The way out is, a future renormalization group analysis. One has to require that 
the solutions do not depend on the cut-off 



dM b (A; m, (3; Q) = 0, or 5 A M b + 5 m M b + 8 p M b = 



(91) 



at a particular scale set by the momentum transfer Q 2 . The renormalization group analysis has thus 
far been done only for asymptotically large Q 2 — > oo. It is a well established fact, even in light-cone 
quantization [0, [l3|, 47, fig] , that asymptotically holds (3(Q) — * (3q / \n{Q 2 / k 2 ) , with a coefficient 



Po well-known from theory and a mass scale k which must be determined from experiment. The 
only difference between the textbooks and the present approach is, that the coupling function (3(Q) 
must be calculable down to momentum-transfer zero. The latter is required when one solves the 
integral equation. Here is the problem, and here we must leave it, since the explicit calculation of 
r(Q,A) and the subsequent renormalization of /3(Q,A) to yield a (3(Q,k) clearly goes beyond the 
scope of the present work. But we have spotted where the problem is. 
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6.4 Retrieving the full wavefunction 



One should emphasize finally that the knowledge of the qq-space eigenfunction ifo is sufficient to 
retrieve all desired Fock-space components of the total wavefunction The key is the upwards 
recursion relation Eq,([30|), i.e. 



n-l 



(n\V) = J2Gn(n\H n \j)(m . (92) 
i=i 

Obviously, one can express the higher Fock-space components (n|^) as functionals of ip by a finite 
series of quadratures (or matrix multiplications). One does not have to solve another eigenvalue 
problem. In order to show this, we ask for the probabilty amplitude to find a gg{2) or a qqg(3) 
state in a particular meson b as an example. With (lj^) = {qq\ipb} the first two relevant equations 
of pi are 



(2|tf) = G 2 (2|fr 2 |l)(l|^) , and (93) 

<3|tf> = G 3 <3|H3|1>(1|*>+G 3 (3|fl-3|2)<2|*) . (94) 

Note that (2\^f) is expressed already in terms of The sector Hamiltonians H n have been dealt 

with in sections ||]and[5]. They can be expressed by chains with the bare interaction V, particularly 
H 2 = T + V + VG 3 V + VG 5 V and H 3 = T + V + VG A V + VGqV + VG 6 VG 5 VG 6 V. With 
the abbreviation Vij = (i\V\j) and the 'selection rules' of Figure § one gets (2|i^2|l) = V23G3V21, 
(3|i? 3 |l> = V31, (3\H 3 \2) = V 32 , and therefore 

<2|tf> = G 2 V2 3 G 3 V 21 (l\y), and (95) 

(3|^) = G 3 y3l(l|*) + ^3^32^2^23^3^21(11*) • (96) 

One should emphasize again that these expressions are exact to all orders in the coupling constant. 
The propagators have to be treated like exposed in section They include the running coupling 
constant and in the solution are independent of lo. Obviously, the procedure is straightforward, 
well-defined, and actually simple. 



7 Summary and Perspectives 

In the present work, one adopts the point of view that most, if not all, properties of a Lagrangian 
gauge field theory are contained in the canonical front-form Hamiltonian, which is calculated in 
the light-cone gauge omitting the zero modes. Periodic boundary conditions allow to construct 
explicitly the Hamiltonian matrix for fixed harmonic resolution. Rows and columns of this matrix 
refer to all possible momentum states of 2, 3, or more particles. All variables are well defined, finite 
and denumerable: The total longitudinal momentum and the harmonic resolution are finite, the 
number of states and the number of sectors are finite, and last but not least every single matrix 
element of the Hamiltonian is finite. One is confronted with the diagonalization of a finite block 
matrix and applies the theory of effective interactions in conjunction with the method of iterated 
resolvents. As shown, one can map the large Hamiltonian matrix without smallness assumption 
onto a matrix problem with much smaller dimensions, to the problem of diagonalizing the effective 
Hamiltonian in the qq-space. It is shown why these two matrices have the same eigenvalues. At 
the end, the periodic boundary conditions are relaxed by the limiting process of going over to the 
continuum limit. 

The diagrammatic representation of the effective interaction as given in Figure looks like 
second order diagrams of perturbation theory. But despite this, they are far from being that: due 
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to the vertex function (which after renormalization is related to the running coupling constant), 
they represent a resummation of pertubative graphs to all orders in the bare coupling constant. 
The effective interaction U is computed explicitly in section |6|. It is our understanding that U + U a 
represents an exact mapping of the full many-body Hamiltonian onto the qq space. Solving for the 
eigenfunctions in this space one can retrieve all other many-body amplitudes of the wavefunction 
in a self-consistent way, as shown in section ^. 

To arrive at this simple result it is crucial to include arbitrarily many gluons. Only then, as 
argued in section || the Russian puppet structure of 'resolvents within resolvents' can be used as 
an argument why the effective interaction is the same in all hierarchies. It is precisely this aspect 
of self-similarity in a gauge theory which ultimately allows for the 'breaking of the hierarchy': In 
the solution, the particles propagate like free particles and all many-body aspects reside in the 
vertex coupling function. One can overstress the point for QED by stating that the generation of 
the simple Coulomb potential requires infinitely many photons. Finally, one should emphasize that 
the present approach needs no 'convergence-improving' external potential as in [19], which once 
introduced is so difficult to get rid-off. The average potential generates itself, in the solution. 

The present work claims to connect a Lagrangian gauge field theory with the eigenvalues and 
eigenfunctions of a bound state equation. We like to emphasize that all steps can be (dis-)verified 
numerically in a well defined manner. One can start with a sufficiently large DLCQ-Hamiltonian 
matrix, as large as the computer can digest. The reduction to the effective Hamiltonian can 
then be done by successive matrix inversions and multiplications as described in Sections |3| and |j. 
The prediction is that the lower eigenvalues depend only weakly on the starting point energy. In 
addition, one can (dis-)verify numerically a posteriori how well the succesive approximations from 
Eqs. ( |57[ ) to ( |58| ) to (|^) are satisfied in practice. 

Important as the present work might be, it only can be the first step. The resulting effective 
potential depends on several formal and unphysical parameters, among them the harmonic reso- 
lution K and the transversal length L± induced both by the periodic boundary conditions, the 
cut-off scale A regulating the transversal momenta, and finally the smallness parameters e and A. 
The resolution and the length disappear when going to the continuum limit. But the dependence 
on A, e and A remains. The second, and perhaps even more important step, must therefore be 
to remove them by a renormalization group analysis [45, [&q| . This has not been done, yet. The 
present work thus culminates and ends at a pre-renormalization stage. The merit of the present 
work is to pinpoint the object of renormalization, namely the vertex function which is given here 
for the first time in a non-perturbative formulation to all orders in the coupling constant. 

But even before a future renormalization group analysis the present formalism is useful and 
has some rather nice aspects: Even without knowing the explicit structure of the running coupling 
constant in the infinite momentum frame one can do some educated guess work. For example, 
using Richardson's version [49] of the running coupling constant one can establish confinement 
in a parameter-free fashion. One can approximately, but analytically, calculate the heavy meson 
masses based on the parametric variation of a trial wave function [ 50 1 . An explicit, but numerical, 
solution of the same problem is currently being attempted fl5l}| . Last but not least one can choose 
the smallness parameters A so small that no transversal momenta survive at all. One then ends up 
with a 'collinear model' [^| similar to the dimensionally reduced models solved thus far [24, 3(| 37], 
resulting in stupendous analytical properties of the light mesons [52]. These features make it worth 
to communicate the present work even prior to a renormalizaton group analysis. 
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A The 4x4 Block Matrix as a Paradigm 



The considerations in section |] can be substantiated with an explicit example. The rows and 
columns of any finite dimensional matrix can be grouped into four blocks. The block matrix has 
then the following shape: 



/(1\H\1) (l\H\2) (l\H\3) (1|#|4)\ 

<2|JT|1> (2|JT|2) (2|fT|3) (2|fT|4) 

(3|iJ|l) (3\H\2) (3\H\3) (3\H\A) 

\ (4\H\l) (A\H\2) (A\H\3) (4|if|4)7 



(97) 



The reduction from the 4 x 4 to the 3x3 matrix is easy: Since = H, one replaces all bare block 
matrix elements in the 3x3 matrix by (i\H^\j), with 

H 3 = H + HdH . (98) 

Next, reduce the matrix of block matrix dimension 3 to the one of block matrix dimension 2, by 
H 2 = H3 + H3G3H3. Inserting H3 from Eq.(|98|), one gets 



H 2 = (H + HG A H) + (H + HG 4 H)G 3 {H + HG A H) . (99) 
Performing reduction and substitution once more one arrives at 

fli = ({H + HG 4 H) + (H + HG 4 H)G 3 (H + HG 4 H)^G 2 

((H + HG 4 H) + (H + HG A H)G 3 (H + HG A H)^j . (100) 

One has thus expressed the matrix elements of the effective interaction in sector 1 in terms of the 
bare interaction H and the resolvents G 2 , G3, and G4, which in turn are given by the effective 
interactions H 2 , H3 and H4, respectively, i.e. 

G 2 = |2) -1— (2\ , G 3 = |3) -J— (3| , G 4 = |4> -J— <4| . (101) 

LO — tl 2 U! — 113 to — tl 

Note that all operations are well defined matrix multiplications and inversions. The bare interac- 
tions H alternate with the resolvents of effective interactions H n to build up strictly finite chains. 
The longest chain in Eq.( |100| ) has 5 propagators: HG^HG3HG 2 HG3HGiH . This is in stark 
contrast to the infinite chains of perturbative series. 



B Continued Fractions in Iterated Resolvents 



The structure of the chains in the Method of Iterated Resolvents depends quite strongly on the 
block matrix structure of the matrices considered, as is obvious from Eq.( |100| ) in Appendix [A|. For 
to show that explicitly, the example of a matrix with a tridiagonal block structure is choosen, i.e. 



H 



B 



H 2 \ 



H\ 2 




■ \ 


H 22 


H 2 3 


H3 2 


H33 


H34 




H43 


Hu I 



or H B 





1 






1 


2 


3 




3 


4 


5 


V. 




5 


6/ 



(102) 



with dots representing zero matrices (Hb) or straight zeros (Hb). Applying the same procedure 
as in Appendix [A|, one notes that many of the long chains vanish since for example G 2 HG^ = 
G 2 H^Gi = 0. The effective interactions, as given in Eq.(|100D for the general 4x4 block matrix, 
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Figure 9: The energy function E(uj) for the 
matrix Hb is plotted versus uj. The solutions of 
E(u>) = to agree with the four exact eigenvalues 
-1.87208, -0,01518, 3.33343, and 10.5538. 



simplify strongly. Here, one gets consecutively H% = H + HG^H, H2 = H + HG3H and H% = 
H + HG2H. The effective Hamiltonian in the 1-space therefore has the structure of a continued 
fraction: 

(l|#i|l) = H n + H X2 1 H 2l . (103) 

u) — H22 — H23 j H32 

& — H33 — i?34 — - i?43 

id — -H44 

One can test the method of iterated resolvents in an almost trivial way: Assume that each of the 
blocks in Eq.( |102| ) has dimension 1. The matrix is thus a usual 4 x 4-matrix. The inverse of a 1 x 1 
'matrix' is trivial. The 'effective Hamiltonian' is also a 1 x 1 'matrix' identical with the energy 
function E(u). For the matrix Hb m Eq. Q102; ) one gets simply 



E(u) = . (104) 

uj-2 



5 • 5 

uj — 



uj — 6 



As displayed in Figure |9|, the solutions of the fixpoint equation E(u>) = uj agree with the four 
eigenvalues to within computer accuracy. Note that a form like Eq.( |103| ) or ( |104|) could possibly 
be useful to diagonalize a tridiagonal matrix of arbitraly large but finite dimension. 

C The Chains of the Effective Interaction 

The example in Appendix [A| suggests that the effective interaction wants to develop finite chains 
of bare Hamiltonian blocks (i\H\j) alternating with resolvents of sector Hamiltonians H n . This 
can be cast into a systematic procedure, as follows. Sequential application of the recursion relation 



Eq.(31) allows to derive the useful relation 



N 

Hn = H + ^ H m G m H m . (105) 

m=n+l 
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Table 1: All chains of the effective interaction in the qq-space Hi with one, two and three propaga- 



tors Gi which are associated with the QCD matrix in Figure \1 



for refer to the three different sums in the defining Eq.(ll 



are enumerated. The square brackets 



# 


H{ > and H\ ' 


# 


Hf\i 3 MM 


# 


H[ 3) [l 3 ,hJ 2 ] 


# 




1 


SG 2 S 


10 


FG 4 VG 3 VG 2 S 


15 


VG 3 SG 5 VG 2 S 


25 


VG 3 VG 2 VG 3 V 


2 


VG 3 V 


11 


FG 6 VG 3 VG 2 S 


16 


VG 3 VG 6 FG 2 S 


26 


FG 4 VG 3 VG 4 F 


3 


FG 4 F 


12 


FG e SG 4 VG 3 V 


17 


FG 4 SG e VG 3 V 


27 


r Lrgr Lr 2 r Cxg-T 


4 


FG 6 F 


13 


FG 6 VG 5 VG 2 S 


18 


FG 4SGqFG 2 S 


28 




5 


VG 3 VG 2 S 


14 


FG % VG 5 SG 3 V 


19 


FG 4 VG 7 FG 3 V 


29 


FGqSG 4 SGqF 


6 


FG A VG 3 V 






20 


FG 6 VG 7 VG 4 F 


30 


FGqVG^VGqF 


7 


FG%FG 2 S 






21 


FG 6 VG 7 FG 3 V 


31 


FG 6 FG 2 VG 3 V 


8 


FG e VG 3 V 






22 


FGq S GqFG 2 s 


32 


FG 6 VG 3 VG 4 F 


9 


FGqSG^F 






23 
24 


FG % VG W FG 3 V 
FGqFGwFG^F 







To prove this, one writes down Hi = H 2 + H 2 G 2 H 2 and H 2 = H 3 + H 3 G 3 H 3 and combines them 
to get Hi = H 3 + H 2 G 2 H 2 + H 3 G 3 H 3 . This is Eq.(|105|) for N=3. The general case is proven 



by induction. Eq.(105) relates the effective interaction in sector n to the bare Hamiltonian H 
and virtual scatterings into the higher sectors m > n. This only upward scattering has important 
consequences for the structure of the chains. How can one classify them? The number of propagators 
in a chain turns out to be a more useful criterium for that than, for example, the order of the 
coupling constant. If a chain has 3 propagators, it will contribute to H^ 3 \ The bare interaction 
has no propagator and thus is the 'chain' H^°> = H. The effective interaction is the sum of all 
possible chains: 

H n = H®+H®+HP)+H<?)+HW+H® + ... . (106) 



The expansion is finite for any finite dimensional block matrix. After its insertion into Eq.( 105 ) 
one reads off upward recursion relations with the general term 

= J2 (H^GtH^ + H^dH^ + ■■■ + H^GiH^ + H^H® ) , (107) 



l>n 

or explicitly for the first four of them 

(108) 
(109) 
(110) 
(111) 

For convenience, the chains up to order 4 are tabulated explicitly: 

= ]T HG h H , (112) 

h>n 





= E 

l>n 




1 




H m 

n 


= E 

l>n 




+ h1 1) g 1 hP 


)• 


tf(3) 


= E 

l>n 


(hPg^P 


+ h\ 1) g 1 h\ 1) 


+ H^ Gi hP), 


tf(4) 


= E 

l>n 


(H^dH® 


+ H^G l H ( i 2) 


+ H < t 2) G l H[ l) + Hf i) G l H < f ) ) 
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= J2 (HG h HG h H + HG h HG h H) , (113) 

h>h>n 

Hn^ = ( HG^HG^HG^H + HGi 1 HGi 2 HGi 3 H+ 



h>h>h>n 



HG l3 HG h HG h H + HG h HG h HG h H j ( 114 ) 



+ J2 ( HG[ 3 HG ni HG r2 H 



r2 >ni >™ 



(115) 



Hn^ — ( HG^HG^HG^HG^H + HGi 1 HGi 2 HGi 3 HGi 4 H+ 

l4>h>h>h>n 

HGi 4 HGi 4 HGi 2 HGi 3 H + HGi 2 HGi 1 HGi 3 HGi 4 H+ 
HGi A HGi 3 HGi x HGi 2 H + HGi 3 HGi 2 HG[ 1 HGi 4 H+ 
HGi 4 HGi 2 HGi 1 HGi 3 H + HGi 3 HGi 4 HGi 2 HGi 4 H 
+ Yl [ HGi 4 HGi 2 HGi a HG ri H + HG ri HGi 3 HGi 2 HGi 4 H 

l 4 >l 3 >l2>n K 
l4>l^>ri>n 

+ J2 ( HGi 3 HGi 4 HG ri HG r2 H + HG r2 HG ri HGi 4 HGi 3 H+ 

l 4 >l 3 >n ^ 
l4>r 2 >ri>n 

HGi 3 HGi 4 HG r2 HG ri H + HG ri HG r . 2 HGi 4 HGi 3 H 

Table EJ enumerates all chains up to order 3 for the full QCD-block-matrix as displayed in Figure El 
For convenience, the type of interaction in the block matrix is accounted for by V, F, and S, refering 
to vertex-, fork-, and seagull-interactions, respectively. 



D The Gauge Trick 

The present work has been started keeping track explicitly of the gauge remnants, the seagull and 
the fork interactions S and F. With the work progressing it was realized that the formal manip- 
ulations can be substantially simplified when one includes only the vertex interactions combined 
with the gauge trick. Rather than formally, the argument is lead here diagrammatically, by way 
of example in Figure |l0] and its caption. Very instructively, one misses there a chain of order 1 
which shows up in Table [l], particularly FGqF. It is found when analyzing the series in Eq.(41) 



term by term, some of which being displayed diagrammatically in Figure y. Consider the graph 
labelled (36b) in the Figure, standing for the chain VG3V Gq VG3V. Now, interprete the two outer 
gluon lines as instantaneous and write down the corresponding chain: it is the missing FGqFI The 
literally hundreds of graphs analyzed diagrammatically in preparing this work, cannot be discussed 
here in detail, of course. Suffice it to state that all chains with instantaneous interactions as enu- 
merated in Table [J have been identified explicitly as the instantaneous partners of dynamical chains 
like VGVGV . . . GVGVGV . This should be sufficient evidence for the 'gauge trick ' applied in this 
work. 

E The cancellation of the gauge remnants 

Substituting the 77-dependent terms of the spinor sum, Eq.(^), into the effective interaction yields 
straightforwardly 



^, -, vpt |'- \/'- ng-l g V(Q,A) (P+)< 

2_^{q,q\VR l W ,q,g)(q ,q,g\RV\q ,q) v = — 5 ft^t ^+ 

^ ^Tlc u L 1 XgrCg 
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1 3 2 3 1 




1 3 2 1 




Figure 10: The annihilation interaction with all instan- 
taneous lines explicitly inserted. The annihilation graph 
of the effective interaction is displayed in the upper left 
of the figure. It represents the following three-step proce- 
dure: (1) A line in between two interaction points is called 
a 'physical' line representing the exchange of a physical 
particle; (2) Every physical line is the sum of a 'dynam- 
ical' and an ' instantaneous' line; (3) Two instantaneous 
lines cannot be connected at a vertex. Drawing all pos- 
sible diagrams with 0,1,2,... instantaneous lines like 
in the left of the hgure in the convention of Tang J§^/, 
one obtains seven gaphs. On the right, the same seven 
diagrams are drawn in the convention of Brodsky p3]] . 
One can read off there the chains of different order, 
particularly SG 2 S, FG 4 F, VG 3 VG 2 S, VG 4 VG 3 F, and 
VG3VG2VG3V . When analyzed in the same way, the 
gluon exchange graph Ux t \ in Figure ^ provides one with 
the chain of order zero (S) and with VG3V. All of them 
appear also in Table ||. But one misses a chain of order 
1 which appears in the Table, i.e. FGqF. Where is it, 
and where are the other chains of higher order? — See 
the discussion in the text. 
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u{q)l^u{q')\ [u{q) lu ri u u{q')] , [u{qh^u(q')} [u{q)luk v g u{q') 



P^ \ I XgXg 



P^~ \ I XgXg 



+ 



P~^~ \ I XqXq 



P^~ \ I XqXq 



(116) 



The well-known property of the Dirac spinors 



{k q -k' q Y u(k q ,\q)i^u(k'\') 



(117) 



can be used for constructing the time-like null vectors derived in section ^, i.e. 



'< "> 2/,-,; 



rf 



— rf 



2kg 



(118) 



Together with the resolvent G3 = —x g /k 2 ± one gets most directly 

q\ Vrf \q', q, g)G 3 (q', q, g\ RV \q' , q') v 



n 2 - 1 g 2 r 2 (Q, A) 1 [u(qh + u(q')] [u(g)7 + -u.(</)] 
2n c nP+ xf _ ^4 ~~~~ 



XqXq \ XqXq 



(119) 



Compare this with the seagull interaction evaluated directly by means of the tables 



(q,q\ rfSR\q',q') 



n 2 c -l g 2 r 2 (Q,A) 



u(k q , X q )-f + u(k f q , X' q ) u(kq, \q)~i + u(k'g, \' g ) 



2n c 



( Xq -x') 2 P+ 



X nX 



q^q 



P^ \ I XqXq 



(120) 



to conclude that all gauge artefacts cancel each other precisely. 
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